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ABSTRACT 



In this paper we present for the first time high spatial resolution K-band maps of the central kinematical and near-infrared spectral 
properties of the giant cD galaxy in the Fornax cluster, NGC 1399. We confirm the presence of a central velocity dispersion dip within 
r < 072 seen in previous long-slit studies. Our velocity dispersion maps give evidence for a non-symmetric structure in this central 
area by showing three cr peaks to the north-east, south-east and west of the galaxy centre. Additionally we measure near-IR line 
strength indices at unprecedented spatial resolution in NGC 1399. The most important features we observe in our 2-dimensional line 
strength maps are drops in Nal and CO (2-0) line strength in the nuclear region of the galaxy, coinciding spatially with the drop in cr. 
The observed line strength and velocity dispersion changes suggest a scenario where the centre of NGC 1399 harbours a dynamically 
cold subsystem with a distinct stellar population. 
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1. Introduction 

The centres of giant early-type galaxies are interesting labora- 
tories for exploring important questions about the physics of 
early-type galaxy assembly and star formation history. In partic- 
ular, the ultimate goal is to tie the formation of global properties 
to such central properties as nuclear stellar populations, super- 
massive black holes, and active galactic nuclei. 

The recently developed technique of near-infrared integral 
field spectroscopy in conjunction with adaptive optics offers a 
new window into these astrophysical questions by probing spa- 
tial scales (~07 1 or better) heretofore only accessible to the in- 
struments on the Hubble Space Telescope (HST). At the distance 
of the Fornax cluster (~ 20 Mpc) this corresponds to a spatial 
resolution of about lOpc, which can be achieved over a 3" x 3" 
field-of-view with e.g., VLT/SINFONI. Although with limited 
field of view, this setup can be used to scrutinise compact stellar 
subsystems such as globular clusters or the nuclear regions of 
galaxies. 

Today it is widely accepted that many massive and 
bright early-type galaxies contain super-massive central black 
holes (SMBH) with masses reaching up to 1O 9 M , which 
tightly c orrelate with the physical p r operties of the host 
galaxies (iKormendy & Richstond [19951 iGebhardt et ail 120001: 
iFerrarese & Merrittll2000l) . Recently, it was shown that the nu- 
clei of low- and intermediate luminosity early-type galaxies 
in the Fornax and Virgo clusters exhibit a (often blue) cen- 
tral (<0.02 effective radii) light excess in contrast with high- 
luminosity galaxies (Mb ^ -20), which usually s how a flat 
luminosity profile (core) at the same spatial scales (Cote et al. 
120061 120071) . The mass fraction contributed by the central light 
excess is similar to the relative contribution of SMBH in gi- 



ant galaxies, which could suggest a common formation path 
(IWehner & Harrisl 120061: fCote et al.ll2006h . The nature of these 
central light anomalies is still uncertain since until recently their 
spatial scales were typically reachable only by HST imaging and 
in a few cases HST/STIS spectroscopy. 

In this paper we present for the first time high spatial resolu- 
tion two-dimensional maps of the velocity profile and the distri- 
bution of two near-IR spectral features in NGC 1399 - the cen- 
tral giant elliptical galaxy in the Fornax clust er (Mb ^ —21.08) . 
This galaxy is known to harbour a SMBH (Sagliaetal. 2000; 
Houghton et alj 20061: Gebhardt et ai1l2007l) . weak nuclear ac- 



tivity (O'Connell et al 



2005) and shows a classical core in the 



* Based on observation collected at the ESO Paranal La Silla 
Observatory, Chile, during the Science Verification of SINFONI. 



luminosity profile dCote et al.ll2007l) . 

The paper is organised as follows: in Section 2 we present 
the observations and data reduction steps. Section 3 is devoted 
to the analysis of the kinematics in NGC 1399, while Section 4 
presents the measurements of near-IR line strength indices. The 
main observed features and potential explanations are discussed 
in Section 5. Our concluding remarks are given in Section 6. 



2. Observations and data reduction 

2.1. Observations 

The observations were made with SINFONI - a near-infrared 
Integral Field Unit (IFU) spectrograph, mounted in the 
Cassegrain focus of Unit Telescope 4 (Yepun) of VLT at Paranal 
La S illa Observatory ( Eisenha uer et al.l l2003t iBonnet et al.l 
2004). The instrument's field of view on the sky is divided into 
32 slices. The pre-slit optics allow us to chose between three 
different spatial samplings: 0725, a/1 or 07025, yielding a total 
field of view of 8" x 8", 3" x 3" or a/8 X a/8, respectively. On 
the detector each slitlet is imaged over 64 pixels in a brick-wall 
pattern and thus each one is well separated from the neighbour- 
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ing slitlets on the sky. The Science Verification detector used 
during our observations had a blemish affecting wavelengths be- 
tween 2.37 and 2.40 fim of the spectrum from slitlet number 17. 
This defect had no impact on our investigation since it lay out- 
side of the wavelength region of interest. 



Table 1. NGC 1399 observing log. 



UT date 


t(s) 


O+S Frames 


Total t (s) 


OB ID 


(1) 


(2) 


(3) 


(4) 


(5) 


2004 Oct 05 


2x300 


4+4 


2400 


A 


2004 Nov 25 


2x300 


4+4 


2400 


B 


2004 Nov 26 


2x300 


4+4 


2400 


C 


2004 Nov 26 


2x300 


4+4 


2400 


D 


2004 Nov 26 


2x300 


4+4 


2400 


E 



Notes: (1) Date of observation, (2) exposure time per object pointing, 
(3) number of Object + Sky frames, (4) total on source exposure time 
per observing block, (5) observing block ID. 



Our study makes use of the 071 scale with a grism cover- 
ing the K-band from 1.95 to 2.45 fim at a dispersion of 2.45 
A/pix. In this mode, SINFONI delivers a spatial sampling of 
0705 x 071. However, during data processing the observations 
were re-sampled to 0705 x 0705 as allowed by our dithering 
strategy. The spectral resolution around the centre of the band 
is R^3800, measured from arc lamp frames. NGC 1399 was ob- 
served on 5 October, 25 November, and 26 November 2004 as 
part of the SINFONI Science Verification program. Observing 
conditions were generally good with free-air optical seeing mea- 
surements ranging from 074 to 077, except in the night of 
October 5 , where they reached 1". 

Observing in the near-IR spectral domain implies several 
challenges. One of them is the very bright night sky consisting 
of many strong emission lines. To remove the night sky signa- 
ture we used the standard near-IR technique of consecutively 
taking object (O) and sky (S) frames. Our target was observed 
in five observing blocks (OBs) each containing a OSSOOSSO 
sequence (see TableQ]). Each individual object pointing had two 
300 s integrations, the sky pointings only one. The various on 
source pointings were dithered by 0705 to reject bad pixels and 
assure better spatial sampling. In the course of the data reduc- 
tion we decided to exclude OB A, because of highly variable 
sky conditions resulting in a poor sky subtraction. 

In order to improve upon natural seeing the observa- 
tions were carried out using adaptive optics (AO) in natu- 
ral guide star assisted mode. We used the s ame guide star 
2MASS J033829 14-3526442) as in the study of lHoughton etail 
2006) which is located 1776 to the North of the galaxy centre. 

Together with the science observations appropriate telluric 
standard stars were observed, as well as four velocity template 
stars covering the spectral range K4III to K7III (see Table|2]). 



2.2. Basic data reduction and telluric correction 

The ob servations were reduce d with the ESO SINFONI Pipeline 
v. 1.6 (Modigl iani et alj|2007l) . Calibration files such as distor- 
tion maps, flat fields and bad pixel maps, were obtained with the 
help of the relevant pipeline tasks ("recipes"). Finally, for each 



Table 2. List of velocity template stars. 



Name 


spectral type 


Source 


(1) 


U) 


(3) 


HD 0005425 


K4III 


this paper 


HD 0025211 


K4/5III 


this paper 


BD-15 1319 


K5III 


this paper 


BD-16 1418 


K7III 


this paper 


2MASS J2041 1845+0016280 


M0III 


075.B-0495(A) 


HD 0141665 


M5III 


075.B-0490(A) 


2MASS J17093801-2718559 


M5II-III 


075.B-0490(A) 



Notes: (1) star name, (2) spectral type, (3) origin of data. 



data set contained in one OB the j itter recipe extracts the raw 
data, applies distortion, flat-field corrections, wavelength cali- 
bration and stores the combined sky-subtracted spectra in a 3- 
dimensional data cube. The same pipeline steps were also used 
to reduce the observations of the telluric and velocity template 
stars. The jitter recipe also provides a one-dimensional spec- 
trum extracted from an optimal aperture, which worked well for 
our star observations. 

The accuracy of the wavelength calibration was checked by 
comparing the position of known sky emission lines to the wave- 
length calibrated sky spectra. Where needed small shifts of up to 
one pixel (= 2.45 A) were applied to our galaxy observations. 

The next step is to correct the observations for telluric ab- 
sorption lines, which are especially deep at the blue end of the 
K-band. These spectral lines may vary with time and position on 
the sky. For each observing block, one solar-type telluric stan- 
dard star at a similar air-mass was observed. Each telluric star 
was divided by a scaled and Gaussian broadened solar spectrum 
to match the resolution of our observations. In this way we re- 
moved the spectral features typical for solar type stars as well as 
ensuring a relative flux calibration. The last step before applying 
the telluric correction was to scale and shift the telluric spectrum 
for each individual data cube by a small amount (< 0.5 pix) to 
minimise the residuals of the telluric lines (for details of the pro- 
cedure see Silva, Kuntschner & Lyubenova 2008). Finally, each 
data cube for the galaxy and velocity template stars was divided 
by the optimised telluric spectrum to yield fully calibrated data. 

In order to achieve the best signal-to-noise ratio we com- 
bined four individual data cubes (OB ID: B, C, D, E). For this 
purpose we identified the centre of the galaxy in each data cube 
with an accuracy of one pixel (= 0.05") and re-centred them. 
The final combination was performed using a sigma-clipping 
pixel reject algorithm. 

The next step was to spatially bin the final data cube to a 
roughly constant signal-to-noise ratio which allows us to mea- 
sure reliable kinematics and line s trengths. This was done us - 
ing the adaptive binning method of Cappellari & Copinl (|2003). 
Using our velocity template stars we first performed Monte 
Carlo simulations to determine the minimum S/N needed to re- 
cover the higher order moments I13 and h.4 of the line-of-sight ve- 
locity distribution (LOS VD) to a precision of about + 0.04. This 
also means that the recession velocity and velocity dispersion 
can be recovered with an error of about ± 17 km s . This accu- 
racy is comparable with or better than what was achieved in pre- 
vious studies, while we also preserve a good spatial resolution. 
For our data set and the spectral region 2.2 - 2.3 /im we found 
that a S/N of about 100 is needed (see also Section |3~TT i. The re- 
sulting bins have a typical size of ~ 0715 (5-6 spaxels combined) 
in the central regions, and approximately 072 (>10 spaxels com- 
bined) for radii >1". Bins which contain more than 30 spaxels (a 



Mariya Lyubenova et al.: Central K-band kinematics and line strength maps of NGC 1399 



3 



few bins at the edges of the field of view) are excluded from the 
analysis, because of higher systematic effects, and are marked 
with black colour on the kinematical and line strength maps. 
During the binning procedure we derive pseudo error spectra by 
determining the standard deviation at each wavelength. While 
this procedure is not fully correct since not all spaxels in a given 
bin are independent (e.g., seeing effects) it takes into account 
other error sources such as sky-subtraction and telluric correc- 
tion residuals. 



2.3. Effective spatial resolution 

A key parameter to characterise our data is the effective spatial 
resolution, which was achieved in the combined data cube be- 
fore binning. Direct measurements of the effective point spread 
function were not obtained during the nights of our observations. 
Therefore, we estimate the effective seeing in two other ways. 

First, we compared the reconstructed and radially averaged 
luminosity profile from our data cube with the results from a 
HST/NICMOS image (GO 7453, PI Tonry; camera NIC2, sam- 
pled at O/0755 per pixel). We found that the two luminosity 
profiles agree very well even without any further broadening of 
the NICMOS data suggesting roughly similar PSF sizes with a 
FWHM =s 0/17 (see Fig.Q}. 

Secondly, we independently confirmed the above estimate 
from our own data cube by measuring the FWHM of a point 
source like object (presumably a globular cluster) located at 
about l'/45 south-east from the NGC 1399 nucleus. Globular 
clusters at the distance of the Fornax cluster might be extended. 
However, here we assume it to be a point source and thus ob- 
tain a conservative estimate for the seeing. A simple Gaussian 
fit yields a FWHM = 0/ 1 1 , very much in agreement with the 
results from the comparison of the luminosity profiles. 

We conclude that the effective spatial resolution of our data 
can be described with having a FWHM of about 0715. 
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Fig. 1. Normalised luminosity profile of NGC 1399 as recon- 
structed from our data cube (solid grey line). As a compari- 
son we show HST/NICMOS F160W (short dashed line) and 
HST/ACS F606W (dotted line) observations (retrieved from the 
HST Archiv e) as well as H S T/WF PC2 F606W (dash-dot line) 
profile from lGebhardt et al l d2007b and a VLT/NACO K-band 
profile from Houghton (private communication, program 078. B- 
0806, PI Houghton) - long dashed line. Profiles are normalised 
in the region 1" < r < VIA. The WFPC2 light profile is super- 
sampled and deconvolved for the instrumental PSF which at 
least partly explains the differences between the WFPC2 and 
ACS data, obtained in the same filter (F606W). 



3. Kinematics 

In this section we present our kinematical analysis of the cen- 
tral region of NGC 1399. Our K-band IFU spectroscopy is used 
to probe the velocity profile of the galaxy's centre to HST like 
spatial resolution. 

3.1. Extraction of the kinematics and comparison with earlier 
studies 

The extraction of the kinematical information from the fi- 
nal galaxy data cube was performed with the help of 
the Penalised Pixel-Fitting method (pPXF) developed by 
Cappellari & Emsellem (2004). The success of this method re- 
lies on the provision of a good set of template spectra, which 
match the galaxy spectra as well as possible. We found that the 
template stars observed in conjunction with the NGC 1399 ob- 
servations (covering K4III to K7III) did not yield fully satisfac- 
tory fits and thus we searched the VLT/SINFONI archive for 
further template stars particularly covering later spectral types. 
We found three suitable stars with spectral type M0III, M5III 
and M5II-III. Fully reduced spectra were kindly provided by M. 
Cappellari (private communication). In our final analysis seven 
velocity template stars (see Table |2| were used in the fitting pro- 
cedure to determine the recession velocity, velocity dispersion 
and the hi, and I14 Gauss-Hermit coefficients for each bin in the 
data cube. 



The fitting was performed in the spectral range between 
2.18//m and 2.35/vm where several prominent absorption fea- 
tures allow fitting (see Fig.|2|. However, we excluded the region 
2.20 //m- 2.22 fim (around the Na I feature), because none of the 
template stars was able to match the strength of the Na I absorp- 
tion in NGC1399. This is not surprising, since lSilva et all {2008) 
showed that the Na I absorption strength is typically stronger in 
early-type galaxies compared with Galactic open cluster stars. In 
the fitting procedure typically only 3-4 stars were used to achieve 
a good fit over the full wavelength extent. The stars which re- 
ceived the highest weighting were of spectral type M0III and 
M5III. 

We derive one-sigma error estimates from the pPXF code, 
based on our error spectra (see Sect. |2.2| i. The validity of our 
error spectra is confirmed since the reduced y 1 values derived 
by the code were typically close to unity. Furthermore, we per- 
formed an empirical error estimation by assuming spherical 
symmetry in the galaxy and deriving the standard deviation for 
v, cr, hj, and h\ measurements in a ring of radius 0/8 to 1". The 
errors derived with this procedure agree well with the pPXF in- 
ternal error estimates and are 17.2km s _1 , 19.4km s , 0.03 and 
0.05 for v, cr, I13 and I14, respectively. 

The kinematics of NGC 1399 have been studied several 
times in the past using long-slit spectroscopy, mainly for the pur- 
pose of determining the mass of the super-massive black hole in 
its centre. Velocity dispersions as high as ~ 450 km s make this 
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Fig. 2. Normalised spectrum of NGC 1399 taken from the cen- 
tral bin of our data. The grey line shows the best fitting com- 
posite template as derived in our kinematics measurements. The 
vertical dashed lines show the extent of the central bandpass of 
the Nal and CO (2-0) indices. The positions of several other ab- 
sorption features are indicated. 



galaxy an extr eme and very interesting case in the cr-black ho le 
mass relation dFerrarese & Merritt 2000; G ebhardt et al.l l2000). 

Sagli a et ail ([2000) used ground based optical observations 
and reported little rotation (< 30km s -1 ) and compatibility 
with the presence of a central black hole of M * 5 x 10 8 M . 
iHoughton etail d2006l) studied NGC 1399 with ground based 
adaptive optics assisted observations in the near-IR using the 
ESO VLT instrument NACO. Their analysis reveals a kine- 
matically decoupled core and double-peaked velocity disper- 
sion across the centre. The modelling requires a central black 
hole with M * 1.2 x 10 9 M Q and a strongly tangentially bi- 
ased orbit dist ri bution . The most recent study is the one by 
iGebhardt et all (l2007h . They analysed HST/WFPC2 optical 
imaging and HST/ STIS spectroscopy in combination with the 
Sagli a et ail J2000) data. They see a dramatic increase in the ve- 
locity dispersion at about 075 on both sides of the galaxy and a 
central drop, which they interpret as tangentially biased orbital 
distribution. Their best fit model requires a black hole with mass 
(5.1 +0.7)xl0 8 M G . 

In Fig. [3] we show our kinematical measurements together 
with data from previous studies. We use radial plots for this 
first-order comparison, since the observations in the earlier stud- 
ies wer e obtained with lon g-slits a t different posi t ion an gles 
(110° i nlSagliaetal.l(l2000l). 1 17° in IGebhardt etafl (120071) and 
5.06° in lHoughton et alJd2006h ). We plot data from the full two- 
dimensional region covered by our observations. 

For radii > Of.' 3 we find in general good agreement be- 
tween our measurements and previously published results. 
H owever, we cannot con firm the large negative h\ values seen 
bv lHoughton et al.1 d2006l) although we also measure negative h\ 
values at these radii. We also cannot confirm the velocity disper- 
sion measurements of Gebh ardt et alj ((2007) for their two outer- 
most data points. 

For radii < 072 our measurements app ear to be between 
the me asurements of Hough ton et alj d2006l) and lGebhardt et alj 
d2007l) . Although the studies do not agree well for the inner ve- 
locity dispersion gradient, all three of them show a pronounced 
drop in the very centre. In this comparison the HST/STIS data 
probab ly have the best spatial resolution (^ 0705; Gebhard t et alj 
(2007)), however, even accounting for seeing differences it is 
difficult to reconcile the data sets. Similarly, we cannot find any 
evi dence in our data for th e extreme, positive h\ values seen in 
the lHoughton et all d2006l) study. 
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Fig. 3. Light distribution and the first four moments of the ve- 
locity profile for NGC 1399. Top plot: Luminosity profile, as 
inferred from our SINFONI data. The solid line represents the 
mean luminosity in 07 1 bins. Bottom panels: radial profile of the 
kinematics. SINFONI observations are marked with solid grey 
circles. Empirical error estimates (for details see Section lXTT l are 
shown with an error-bar and a solid symbol in the top-right cor- 
ner of each plot. The solid line represents the mean value of each 
moment in 071 bins. Results from earlier studies are shown with 
different sy mbols as explained on the top plot. Representative 
errors from IHoughton et all (2006) are shown with an error bar 
in the top-ri ght corner o f each pane l with the triangle sy mbol. 
Errors from lSaglia etaTI (120001) and IGebhardt et all d2007b are 
plotted with the data points. 



3.2. Kinematical maps 

Given that we have full two-dimensional kinematic maps, we 
can explore if there is a more complex kinematical structure in 
the centre of NGC 1399, which might explain the measurement 
differences reported above. In Fig. [4] we present our maps of cr, 
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v, fly and hi,. We also ov erpl ot the slit orientations and widths 
of iHoughton et alj d2006l) andlGebhardt et ail (120071) studies for 
better illustration. The lSaglia et al.l (l2000h slit covers almost our 
full field of view, thus it is not shown in Fig. [4] 

The velocity dispersion map shows rising values towards 
the centre with a distinct drop in the middle. However, the in- 
ner ~(X'5 appear to be asymmetric, and show three peaks to 
the nort h-east, south-east and west of the centre. Fig. [4] shows 
that the iGebhardt et all d2007l) slit intersects two of the high-<x 
features, which can p otentially explain the disagreement with 
Houghton et alj (120061) . whose slit does not cover any of them. 
If confirmed, the velocity dispersion structures presented in this 
work may imply the need for more detailed dynamical modelling 
for the nuclear parts of the galaxy. 

Our velocity map does not show any significant ordered mo- 
tion within our errors (^17kms _I ). The h$ map is consistent 
with a value of zero and thus no structure, while the h\ map is 
mildly positive in the very centre and mostly negative for radii 

>a/4. 



4. Near-IR line strengths 

In this section we discuss the two-dimensional distribution of 
two near-IR absorption features - Nal close to 2.2 jum and 
CO (2-0) at 2.3 fim (see Fig. [2j - at a spatial resolution and scale 
only accessible by HST observations in the past. 

4. 1. Index definitions and measurements 

Index defin i tions for Nal and CO (2-0) were taken from 
Fro g el etail d200ll) . The line strength of Nal is measured in 
A, identical to th e method used by t he optical indices in the 
Lick/IDS system dWorfhev et al.l lT994). where the central band- 
pass is flanked to the blue and red by pseudo-continuum band 
passes. CO (2-0) is measured in a similar way but due to the lack 
of a well defined c ontinuum on the red side of the CO feature, 
Fro g el etail ((2001) define four pseudo-continuum band passes 
on the blue side. Thus the continuum used to derive the value of 
the CO (2-0) index is an extrapolation from the blue wavelength 
region. 

The observed spectrum of a galaxy is the convolution of the 
integrated spectrum of its stellar population(s) by the instrumen- 
tal broadening and line-of-sight velocity dispersion (LOSVD) 
of the stars. These effects broaden the spectral feature and in 
general reduce the observed line strength compared to the in- 
trinsic value. In order to compare our observations with other 
studies, often obtained with different instruments, and eventu- 
ally with population synthesis models, one needs to calibrate 
the data to a common system. Since an averaged spectrum of 
the observations reported in this study was already used in the 
ISilva et al.l d2008l) VLT/ISAAC study of early-type galaxies, we 
decided to use that spectral resolution here as well. Therefore, all 
galaxy spectra were broadened to a spectral resolution of 6.9 A 
(FWHM) matching the resolution of ISAAC in the K-band. 

Using our measurements of the velocity for each bin (see 
Sect. I3.lt we can determine line strength values. However, be- 
fore we can analyse the line strength maps we need to account 
for line broadening due to the internal motion of stars in the 
galaxy. In order to compare the index measurements for regions 
with different LOSVD or between different galaxies, we cali- 
brate the indices to a zero velocity dispersion. Our LOSVD cor- 
rections were derived by broadening the velocity template stars 
to velocity dispersions ranging up to 440 km s and h^, ha, val- 



ues between -0.20 and +0.20. For each index we pa rametrised 
the LO SVD corrections following the method by iKuntschnerl 
(|2004|) . Using our own kinematical measurements we calculated 
for each bin and index a correction factor which is applied to 
the raw line strength measurements. At cr — 400km s (h^, 
h\ - 0.0) the correction factor is 1.53 and 1.16 for Nal and 
CO (2-0), respectively. The typical error on the correction fac- 
tor as derived from the differences between the template stars is 
3.3% for Nal and 1.3% for CO (2-0). 

Using our error spectra we compute index errors via Monte 
Carlo simulations, where we take into account photon noise, re- 
cession velocity errors, and LOSVD correction errors. Following 
the procedure applied for the kinematics (see Sect. I3.lt we also 
derive empirical error estimates by computing the standard de- 
viation of index values in a region with radius between / /8 
and 1". For the Nal index, the empirical estimate is 0.5 A, and 
thus slightly higher than our internal estimate of 0.38 A. For the 
CO (2-0) index we find a much larger difference with an empir- 
ical error of 0.6 A as compared to 0.17 A for the internal value. 
Line strength indices can be very sensitive to continuum shape 
effects, especially if they cover large wavelength ranges as is the 
case for the CO (2-0) index. We consider th is the main reason 
for increased errors (see also lSilva et al.l2.QQ8b- In the following 
we adopt the empiric al error estimates in our analysis. 

ISilva et ail d2008l) present a more complete set of spectral 
features in the K-band, including the Ca I and two Fe indices (see 
Fig. |2j. Here we will not discuss these additional, weak features, 
because our chosen SNR per bin in combination with the very 
high velocity broadening of NGC 1 399 makes them very difficult 
to measure. 

4.2. Line strength maps 

In Fig. [5] we show the distribution of Nal and CO (2-0) over 
the central 3" x 3" of NGC 1399. In the central 072 we see a 
well pronounced drop of about 1 A and 1 .5 A in Na I and CO (2- 
0), respectively. Below the maps we show the radial distribu- 
tion of both indices for better illustration. The red line shows 
the median value of ea ch index in 071 wide bins. We note, that 
Hou ghton et alj d2006l) also found some evidence for a central 
CO drop in their data, although their measurements were not 
conclusive. Outside the nuclear region we find evidence for line 
strength gradient in Na I. The index decreases with radius with 
a slope of -0.498 + 0.087 within 074 and 174 (the Spearman 
correlation coefficient is -0.21 with a probability of 0.09% that 
there is no correlation). The linear fit to the data is overplotted 
in Fig. [5] For CO (2-0) we do not find a significant line strength 
gradient over the observed radial range. 

5. Discussion of the central velocity dispersion and 
line strength drops 

In this section we will discuss the significance of the central cr 
and line strength drops we described in Section l3~2l and I4.2l and 
explore a few possible scenarios that could explain their origin. 

5.1. Technical aspects 

Firstly, we investigate the statistical significance of the cen- 
tral drops in absorption strength of the Na I and CO (2-0) in- 
dices. Using a two-dimensional Kolmogorov-Smirnov test we 
can reject the null-hypothesis that the innermost five data points 
(r < 0715) are drawn from the same distribution as the data 
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Fig. 4. Maps of the first four moments of the velocity profile in the central 3" of NGC 1399. Top left - veloci ty dispersion, top righ t 
- velocity. Bottom left - and bottom right - ^4. North is up, east - t o the left. Slit orientatio n and width of the lGebhardt et al.1 (|2007) 
study are shown with solid lines (PA =117°, width = 072), the slit of lHoughton et all d2006l) is shown with dashed lines (PA = 5.06°, 
width = 0/172). Overplotted are levels with constant surface brightness, as derived from our data. 



points between 073 and 1" with a significance level of 3.4% and 
0.4% for the Na I and CO (2-0) indices, respectively. While for 
the Na I index the significance is only marginal, both index drops 
taken together are believable. Another test we performed is to fit 
a simple two component model, where we fit the data with two 
linear functions with a break radius at r = 0/25. The steep slopes 
measured for the central component are different from zero at 
four and nine sigma level for the Na I and CO (2-0) indices, re- 
spectively. 

Secondly, we explore the possibility that the line strength 
drops could be caused by wrong kinematics (e.g., too low ve- 
locity dispersion) which would lead to wrong (e.g., too low) 
LOSVD corrections. Due to the large velocity broadening seen 
in the centre of NGC 1399 we have to apply significant LOSVD 
corrections (see Sect. @J. However, even if we assume the ex- 
treme case of a LOSVD (cr = 370km &-\h 3 = 0.0 and h 4 =0.0) 
constant with radius we still find highly significant drops for both 
indices. 

In summary we find that the line strength drops are signifi- 
cant at least at the 3-<x level and are unlikely to be produced by 
technical/calibration problems or by a chance statistical effect. 
Furthermore the observed features in the velocity dispersion and 
line strengths are present also when one performs the same anal- 
ysis with the individual data cubes before combination, although 
the noise level raises. 



5.2. Astrophysical explanations 

Central velocity dispersion drops hav e so far been ob- 
serve d mainly in spiral g alaxies (e .g. iBottema & Gerritsenl 



1997 



Emsell em et a I U iMarauez et all [2003; 
Falcon-Barroso et al . 2006). cr-drops h ave also been found 

(e.g. iGraham et al.1 fl 998; 
20031: lEmsellem et al.l 



in a few early 



type galaxies 
2002; Pinknevetal. 



Simien & PrugnieJ 
2004). For spiral galaxies iWozniak et al.l Q2003) suggested that 
these cr-drops are the result of gas accretion followed by star 
formation and are therefore indicative of a recent in-fall of dis- 
sipativ e material. In a follow up study Wo zniak & Ch ampavert 
(2006) showed that the life time of such cr-drops can be long 
(> 1 Gyr), however this requires continuous star formation in 
the nuclear region at least on a level of 1 M Q yr~ l . Without con- 
tinued star-formation, the amplitude of the cr-drops decreases 
rapidly. 

Velocity dispersion drops were found in about half of the 
sample of spiral galaxies, studied by Peletie ret al.l d2007l) . Their 
stellar populations analysis, however, shows that the ages of 
these inner components are not significantly younger than the 
centres of spiral galaxies without cr-drops as would be expected 
from the above scenario. This may indicate that the central, dy- 
namically cold components can be more long lived. 
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Fig. 5. Top plots: Maps of Na I and CO (2-0) in the central 3"of NGC 1399. Overplotted are levels with constant surface brightness, 
as derived from our data. Bottom plots: Radial distribution of the same indices. The red lines represent the median index values in 
07 1 bins. The black line shows a linear fit to the index values of Na I in the region 074 < r < 174. The slope value is given in the 
plot. Typical index errors (see Sect. 14. 11 ) are shown with error bars in the upper right corner of the plots. 



Typically the observed cr-drops in spiral galaxies have spa- 
tial scales of a few hundreds p arsecs. For NGC 13 99 such a large 
scale cr-drop is not observed dGraham et al.ll 1 998h . The cr-drop 
we observe in NGC 1399 has a diameter of only -30 pc (cor- 
responding to 073). Similar size cr-drops have been seen with 
HST/STIS spe ctroscopy in the ellip tical galaxies NGC 4649 and 
NGC 4697 by iPinknev et all d2003l) indicating that these small 
scale drops may be a more wide spread phenomenon in early- 
type galaxies. 

On the basis of their dynamical modelling Geb hardt et alj 
(2007) state that NGC 1399 possesses extremely tangentially bi- 
ased orbital distribution for radii between 071 and 075 which 
is responsible for the cr-drop, consistent with the findings of 
iHoughton et alj d2006l) . Both studies explore various possibili- 
ties for the progenitor of these extreme orbits, ranging from nu- 
clear activity, an eccentric disc, a torus to a stellar cluster having 
fallen in on a purely radial orbit. 

We can empirically test the hypothesis for the presence of a 
low velocity dispersion subsystem in the nucleus of NGC 1399 
causing the cr-drop. We performed simple simulations by com- 
bining an averaged galaxy spectrum taken at radii between 072 
and 073 (cr ~ 400 km s ) with a velocity template star, broad- 
ened to velocity dispersions from 15 to 150km s and a range 
of luminosity weightings. We can reproduce our measurements 
for the inner 0715 of NGC 1399 (cr < 360 km s" 1 and positive 
hi\) with ~15% of the light originating from a central stellar sub- 
system with cr up to 100 km s . Using the surface brightness 
profile of NGC 1399 (e.g. from lGebhardt et al.|[2007h we can es- 



timate the surface brightness of such a subsystem. We find that a 
central, low velocity dispersion (cr < 100 km s~') stellar subsys- 
tem with V - 20.3'" is able to explain the central drop in velocity 
dispersion for NGC 1399. Such a V-band luminosity is compa- 
rable to the brightest globular clusters in NGC 139 9 or bright 
nuclei of dwarf galaxies in the Fornax cluster (e.g. lDirsch et al.l 
120031: lMieske"etani2006l) . 

If such a subsystem is present in the nucleus, we can expect 
that it will not have the same stellar population properties as the 
underlying galaxy since it might have been formed in different 
time scales, under different physical conditions or in a different 
environment. Thus we can expect that it will have some influ- 
ence on the line strengths in the same region. Indeed, as shown 
in Sect. 14.21 we find line strength changes with the same extent 
as the cr-drop. Nevertheless, index changes might be due to dif- 
ferent reasons than stellar population effects. A plausible expla- 
nation for the reduced line strengths of Na I and CO (2-0) in the 
central regions could be the prese nce of nuclear activ ity causing 
the continuum level to be raised. lOliva et ail {T999) show that 
the nuclear continuum in the near-IR spectral domain in galax- 
ies with active galactic nuclei (AGN) is mainly due to repro- 
cessed radiation, coming from the AGN UV-continuum heated 
dust. This radiation has a black body spectrum with tempera- 
tures typically between 100 and 1 000 K. |Q' Connell et all ([2005 ) 
report an ultraviolet light outburst in the nucleus of NGC 1399 
with a maximal far-UV luminosity of 1.2 x 10 39 erg s _1 reached 
in January 1999. Thus they claim NGC 1399 to harbour a low- 
luminosity active galactic nucleus (LLAGN). We tested this hy- 
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pothesis by subtracting black body spectra with different temper- 
atures in the above mentioned range from our innermost spectra 
with the goal to match the line strengths seen at radii > 072. 
The results from this test were unfavourable for a scenario in 
which the LLAGN induced non-stellar continuum causes the ob- 
served line strength changes, since we were unable to model 
the index changes for both Na I and CO (2-0) simultaneously 
with one single black body spectrum of a given temperature. 
We can go even further and force the indices to drop by the 
right amount by subtracting a continuum with a very steep slope 
(in comparison with the already mentioned black body spectra). 
However, when we use this steeper slope to match the index 
drops, the effective velocity dispersion measured from such a 
spectrum (cr = 384km s -1 ) does not agree well with our results 
from Sect. 13. II Moreover, the resulting spectrum does not resem- 
ble the observed central galaxy spectrum (i.e. different contin- 
uum shape and slope), which it is supposed to model. Therefore, 
we conclude that the presence of an AGN-like continuum is not 
very likely to be the reason for the central line strength and ve- 
locity dispersion changes observed in NGC 1399. 

Coming back to the s tellar population issue it is very in- 
teresting to note that the Gebha rdt et alj d2007|) imaging data 
of NGC 1399 (HST/WFPC2 PC chip) show a blue nucleus 
(r < 071) with weak evidence for a central light excess at 
the same radii (see FigQ]). Their (B- V) colour drops by about 
0. 1 mag in comparison with the surrounding region, which ap- 
pears to be constant in colour (their Figure 1, bottom panel). 
Recent HST/ACS observations of NGC 1399 in g and z bands 
(Cot e et al.ll2007l also Andres Jordan, private communication) 
show very weak evidence for a central light excess and also 
very mild evidence for a bluer colour, but exhibit prominent 
PA and ellipticity changes within r < 072 consistent with the 
HST/WFPC2 data. We note, that the core region of NGC 1399 
with a break radius of about 2" (see Fig. |T} is much larger than 
the scales discussed here. 

The information from the imaging, combined with the ob- 
served cr and line strength drops at similar radii provide good 
evidence for the scenario of a central dynamically cold stellar 
subsystem with different stellar populations as compared to the 
surrounding galaxy. 

Until stellar population models in the near-IR become avail- 
able we can only perform a limited stellar p opulation analysis 
by using the empirical relations found in the Si lva et all |2008) 
study (their Figure 14). For early-type galaxies with old stel- 
lar populations they found a good linear correlation between 
the Nal index and the optical metallicity indicator [MgFe]' as 
well as CO (2-0) and [MgFe]'. If the central index drops are 
caused by a metallicity change alone (i.e. lower metallicity), 
we would expect to s ee near-IR i ndex c hanges consistent with 
the relations seen in ISilva et all d2008l) . However, this is not 
the case. The index change seen for Nal is much too weak 
compared to CO (2-0), or vice versa. ISilva et alj d2008l) show 
also that centres of Fornax ellipticals with optical signatures 
of a younger population (i.e. strong H/3, indicating luminosity 
weighted ages of around 2Gyr) have a stronger Nal index in 
comparison with old galaxies at the same central velocity dis- 
persions. Also CO (2-0) is expected to increase for intermediate 
ages (~ 1-2 Gyr) since the fraction of thermally pulsating asymp- 
totic giant branch star s near or above the tip of the first ascent 
giant branch goes up (Mara stonll2005l and references therein). 
A scenario w hich is then more consistent with the empirical 
relations from I Silva et all d2008l) is a drop in metallicity com- 
bined with a younger age where the Na I index reacts stronger to 
the age change than the CO (2-0) index. Alternatively, one can 



consider - in addition to a drop in metallicity, abundance ratio 
differences between the main body of NGC 1399 and the pro- 
posed central stellar subsystem, which would explain the relative 
changes of CO (2-0) and Na I. While a combination of metallic- 
ity, age and abundance ratio changes are an attractive explana- 
tion for the observed line strength maps, a confirmation of the 
above scenarios clearly awaits the more detailed predictions of 
stellar population models. 

6. Concluding remarks 

In this paper we present for first time high spatial resolution 
K-band maps for the kinematical and near-IR spectral proper- 
ties of the centre of the giant cD galaxy in the Fornax cluster, 
NGC 1399. We confirm the presence of a central velocity disper- 
si on dip within r < 072 prev iously seen in th e long-slit studies 
of lHoughton et all d2006l) and lGebhardt et all d2007l) . Our veloc- 
ity dispersion map gives evidence for a non-symmetric struc- 
ture towards the centre by showing three cr peaks to the north- 
east, south-east and west of the galaxy centre. Such a complex 
structure in the velocity dispersion has not been seen earlier in 
other early-type galaxies. The west - south-east peaks seem to 
be symmetric about the galaxy centre, while the north-east peak 
is essentially off-centred, nevertheless within our error analy- 
sis we consider it to be real. We note that the nuclear region of 
NGC 1399 does not seem to be obscured by dust (Andres Jordan, 
private communication). 

Additionally we measure two near-IR line strength indices 
at unprecedented spatial resolution. The most important features 
we observe in our 2-dimensional line strength maps are drops 
in Na I and CO (2-0) line strength in the nuclear region of the 
galaxy, coinciding spatially with the drop in cr. 

As an explanation for the line strength and velocity disper- 
sion drops observed in the central region we discuss a scenario 
where NGC 1399 harbours a dynamically cold subsystem with 
a distinct stellar population. Dynamically this can be realised by 
a central stellar disc or a globular cluster (even a dwarf galaxy) 
havin g fallen into the cen tre of NGC 1399 on a purely radial or- 
bit dGebhardt et al. 2007). The near-IR index changes observed 
in our data suggest that the spectral properties of this subsystem 
are directed by some mixture of metallicity, age and abundance 
ratios effects. More precise constraints on the nature of the stel- 
lar population within this cold component awaits more detailed 
stellar population models. 
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